
function DSd=DB_vg_NOlearn(Z,ZZ,v,ep_sd,n)

N=length(v);

R=exp(-ZZ);
V=diag(v.*ep_sd);
Q=V*R*V;
P0I=kron(eye(2),tril(Q)+tril(Q,-1)');

if n>0
    DV=diag([zeros(n-1,1);ep_sd(n,1);zeros(N-n,1)]);
    DP0I=kron(eye(2),V*R*DV+(V*R*DV)');
else
    DR=-Z.*R;
    DP0I=kron(eye(2),V*DR*V);
end

DSd=0.5*diag(P0I).^(-1/2).*diag(DP0I);

